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Heavy quarkonia production is predicted to be sensitive to the formation of the quark 
gluon plasma in relativistic heavy ion collisions via competing mechanisms such as color 
screening and/or quark recombination. During 2004 and 2005 RHIC data taking periods, 
the PHENIX collaboration has measured J/ijj decay into lepton pairs both at mid and 
forward rapidity in Au+Au and Cu+Cu collisions at ^snn = 200 GeV. We present the 
J/ip nuclear modification factor and the J/ip mean square transverse momentum as a 
function of the collision centrality for both systems, as well as the rapidity dependence 
of the J/ip yield for different centrality classes. It is compared to different theoretical 
predictions. All Au+Au and Cu+Cu results shown here are preliminary. 

I. Introduction 

It has been established at lower energies that the J/ips created in nuclear collisions are 
suppressed when the system size and the centrality of the collision increase [|Tj. For heavy 
enough nuclei the suppression exceeds the normal nuclear absorption observed with light 
nuclei. Several mechanisms, including color screening in the quark gluon plasma, have 
been proposed to explain this abnormal suppression observed at low energy. Measurements 
performed at higher energy are therefore crucial to constrain the models and discriminate 
between them. Moreover they may uncover new mechanisms such as quark recombination. 



2. Experimental setup and datasets 

The PHENIX experiment measures J/ip — > e + e~ decay at mid rapidity < 0.35) 
and — > decay at forward rapidity (|?7| G [1.2,2.2]). Electrons are identified 

using RICH detectors and by matching their momentum measured in drift chambers with 
the energy deposited in electromagnetic calorimeters. Muons are selected using a thick 
absorber located close to the interaction point, tracked using cathode strip chambers and 
triggered on using a succession of Iarocci tube planes and steel walls. 

The results presented here correspond to 240 jib^ 1 Au+Au collisions collected in 2004, 
and 3.1 nb _1 Cu+Cu collisions in 2005, both at y^s/vw = 200 GeV. For each centrality, 
rapidity or transverse momentum bin, the J/ip yield is obtained by counting the number of 
unlike-sign dilepton pairs in a mass window centered on the J/ip mass after the background 
has been subtracted using either an event mixing technique or the mass distribution of the 
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like-sign pairs. It is corrected by the efficiency and acceptance of the spectrometer as well 
as the trigger efficiency and normalized to the number of recorded collisions. This yield 
enters the numerator of the nuclear modification factor which is the ratio between the 
measured yield in A+B collisions over the expected yield when assuming binary scaling 
from p+p collisions. The J/ip yield measured in p+p collisions is taken from PHENIX 
2003 published results Systematic errors have been assigned to each contribution. For 
Au+Au collisions, the dominant systematic error comes from the background subtraction, 
due to the poor signal over background ratio, especially for central collisions. 



3. Results 

Figures and 121 represent the J/ip nuclear modification factor as a function of the 
number of participant for Au+Au and Cu+Cu collisions both at mid an forward rapidity 
together with the published measurement from d+Au [ 2 a . A suppression of about a factor 
3 is observed for the most central collisions. 

On figure n the data are compared to R. Vogt prediction assuming 3 mb nuclear 
absorption on top of EKS98 shadowing. The prediction miss the most central points. 
Moreover, from [ 2J it appears that 3 mb is an upper bound for the nuclear absorption 
cross-section in d+Au collisions. 



R. Vogt, nucl-th/0507027, y=2 
R. Vogt, nucl-th/0507027, y=0 
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Figure 1. J/ip nuclear modification factor as a function of the number of participant, 
compared to R. Vogt predictions for the normal nuclear absorption. Vertical bars are 
statistical errors, brackets are point to point systematics and boxes are global systematics. 



On figure El the same data are compared to two sets of models. They all involve 
additional final state interactions and reproduce both the CERN SPS results [ and 
the PHENIX low statistics results from 2002 [0]. The predictions plotted on the left 
panel ([EIEIIZI) overestimate the J/ip suppression. The predictions plotted on the right 
panel show better agreement with the data or even underestimate the suppression. The 
predictions on the right involve either quark recombination mechanisms [ |HJ El E] or 
detailed J/ip transport in medium [CHI- Note however that these predictions differ in 
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the way the cold nuclear absorption is accounted for, in the p+p J/ip production cross- 
section used for normalization and in the open charm production cross-section entering 
recombination mechanism. 
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Figure 2. J/ip nuclear modification factor as a function of the number of participant 
compared to various models of final state interaction in the medium. 



The figure El represents the J/i/j invariant yield BdN/dy as a function of the J/ip rapidity 
for different centrality bins in Au+Au, Cu+Cu and p+p collisions. The vertical lines are 
statistical errors, the bands are point-to-point systematics. Within the large error bars, 
no significant change in the shape of the distribution is observed from p+p collisions to 
most central Au+Au. 
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Figure 3. J/ip yield as a function of rapidity and centrality. 

The figure |U represents the J/i[> mean square transverse momentum as a function of the 
number of collisions for Au+Au, Cu+Cu, d+Au and p+p collisions. Thews predictions [ 
11 without (dashed lines) and with recombination (solid lines) are also shown. 
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Au+Au |y|< 0.35 
Cu+Cu |y|< 0.35 
d+Au |y|< 0.35 
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Figure 4. J/ip mean square transverse momentum as a function of the number of collisions, 
at mid rapidity (left panel) and forward rapidity (right panel). 



4. Summary 

PHENIX measured J/ip in Au+Au and Cu+Cu collisions at y / J/v r /v = 200 GeV at both 
mid and forward rapidity. The dependence of the J/ip nuclear modification factor on 
the number of participants allows to discriminate between some of the models capable 
of reproducing the suppression seen at lower energy. The J/ip rapidity spectra show no 
obvious variation in shape. Interpretation of the J/ip mean square transverse momentum 
is unclear due to large errors. 
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